Role of conviction in non-equilibrium models of opinion formation 
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We analyze the critical behavior of a class of discrete opinion models in the presence of disorder. 
Within this class, each agent opinion takes a discrete value (±1 or 0) and its time evolution is ruled 
by two terms, one representing agent-agent interactions and the other the degree of conviction or 
persuasion (a self- interaction). The mean- field limit, where each agent can interact evenly with any 
other, is considered. Disorder is introduced in the strength of both interactions, with either quenched 
or annealed random variables. With probability p (1-p), a pairwise interaction reflects a negative 
(positive) coupling, while the degree of conviction also follows a binary probability distribution (two 
different discrete probability distributions are considered). Numerical simulations show that a non- 
equilibrium continuous phase transition, from a disordered state to a state with a prevailing opinion, 
occurs at a critical point p c that depends on the distribution of the convictions, the transition being 
spoiled in some cases. We also show how the critical line, for each model, is affected by the update 
scheme (either parallel or sequential) as well as by the kind of disorder (either quenched or annealed) . 
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I. INTRODUCTION 

In the last decades, diverse questions of social dynam- 
ics have been studied through statistical physics tech- 
niques. In fact, simple models allow to simulate and un- 
derstand real problems such as elections, spread of infor- 
mation, vehicle traffic or pedestrian evacuation, amongst 
many others 0. As a feedback, these issues are attractive 
to physicists because of the occurrence of order-disorder 
transitions, scaling and universality, among other typical 
features of physical systems. 

Concerning the particular subject of opinion dynam- 
ics, several models have been proposed in order to study 
the emergence of consensus (for a recent review, see [l|). 
As concrete examples, let us mention opinion models 
based on outflow dynamics majority rules [H-Il] and 

bounded confidence Q, as well as kinetic exchange PT 




1131 ] . Recently, the effects of negative interactions ^ 
and network dynamics [l5l - [T7j in opinion formation have 
also been considered. 

In this work we introduce the degree of persuasion or 
conviction of the agents. It is mimicked by a parameter 
which gauges the tendency of an agent to hold its opin- 
ion or (if negative) change mind spontaneously. Then, 
we study the impact of persuasion in the critical behav- 
ior of a non-equilibrium model of opinion formation with 
a finite fraction of random negative agent-agent interac- 
tions. We study two classes of disorder (either quenched 
or annealed) both for the strength of convictions and for 
agent-agent couplings. We also consider two different 
kinds of update, either sequential or parallel. Numerical 
Monte Carlo simulations show that a continuous order- 



disorder phase transition, where order is characterized by 
a dominating opinion, can occur in all the variants of the 
model considered. However, the critical line is strongly 
affected by the distribution of convictions. Moreover, it 
is also affected both by the update scheme and by the na- 
ture of the random variables, as occurs in other models 

dm. 

This work is organized as follows. In Section II we 
present the opinion model and define its microscopic 
rules. Numerical results are discussed in Section III in 
connection with the analytical considerations presented 
in the Appendix. Section IV contains the conclusions 
and final remarks. 



II. THE MODEL 

We consider an opinion model based on kinetic ex- 
change fToMTs! . At a given time step t, each agent i 
has a discrete opinion Oi(i) = —1,0 or 1, that evolves 
according to 



Oi(t + 1) = Q Oi(t) + fiij oj(t) , 

Oj(t + 1) = Cj Oj(t) + flji Qj(t) , 



(1) 
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where Ci is the conviction of agent i and fj,{j is the 
strength of the influence it suffers from a randomly cho- 
sen agent j in a fully-connected graph. If the value of 
the opinion exceeds (falls below) the value 1 (—1), then 
it adopts the extreme value 1 (—1). Pairwise interaction 
strengths are random variables distributed according to 
the binary probability density function (PDF) 

F(iHj)=p5(jHj + 1) + (1 ~p)S( tHj - 1) . (2) 

In other words, the agents can exchange opinions with 
positive (+1) or negative (—1) influences, and p quanti- 
fies the mean fraction of negative ones [111 ] . In magnetic 



2 



systems, analogous positive (negative) interactions would 
correspond to ferro (anti-ferro) couplings. The influence 
of an individual over another one needs not be reciprocal 
(i.e., not necessarily faj — fj,ji), however, whether inter- 
actions are symmetric or not, does not affect the results. 
If Ci = 1 for all i (i.e., q — 1), one recovers the model of 
Ref. [Ill ], for which there is an order-disorder transition 
at a critical value p c = 1/4. As discussed in Ref. [HI, the 
effect of negative interactions is similar to that produced 
by Galam's contrarians in opinion models 24j . 

However, more realistically, the degree of conviction 
needs not be unitary nor homogeneous |'.KS . Then we 
considered two discrete alternatives for the PDFs of the 
convictions Ci, namely, 
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(3) 
(4) 



They model the cases where a mean fraction 1 — q of 
the individuals have either no convictions or completely 
change mind, respectively. In both cases, the model of 
Ref. [1J| is recovered for q = 1. Notice that, differently 
from the Sznajd dynamics [2J, where each agent inter- 
acts with a group of individuals at a time, in the present 
exchange model, interactions are pairwise. 

We will show that the heterogeneity of convictions may 
provoke the destruction of the order-disorder phase tran- 
sition. Moreover, we will analyze two distinct kinds of the 
random variables Ci and fiij : they can be cither quenched 
or annealed. The former are drawn from the PDFs given 
by Eqs. © and © [or © and (0J] at the beginning of 
each simulation and remain fixed during the evolution of 
the system, whereas the later are renewed at each Monte 
Carlo step (MCs), where one MCs corresponds to N it- 
erations of Eq. ((T|), being N the population size. 

In addition, we have studied two kinds of upgrades: 
synchronous (parallel) and asynchronous (random se- 
quential). In the former case, we randomly choose N 
pairs of agents that interact by means of Eq. ([1]). Only 
after the N interactions took place, the states of the N 
agents are simultaneously renewed, increasing time by 
one MCs. In the asynchronous case, also N pairs of 
agents that interact by means of Eq. |T]) are randomly 
chosen at each MCs, but the opinions are assigned a 
new value at each interaction. A more realistic dynamics 
probably proceeds in between both schemes. 

All simulations start with a random initial distribution 
of opinions, and all interacting pairs of agents are ran- 
domly chosen among the N individuals in the population 
(which corresponds to a mean-field approach). 




(5) 



where (...) denotes disorder or configurational average. 
Notice that O plays the role of the "magnetization per 
spin" in magnetic systems. In addition, we also consider 
the fluctuations \ of the order parameter (or "suscepti- 
bility") 



X = N((0 2 )-(Of) 
and the Binder cumulant U, defined as 

(O 4 ) 



U = l- 



(6) 



(7) 



In the following subsections, we will analyze separately 
the distributions given by Eqs. ^ and (0]). 



A. Model with distribution G\ 

For the distribution Gi(Cj) of Eq. (J3J), the mean frac- 
tion of null convictions is 1 — q. Such agents with no 
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III. RESULTS 

We analyze the critical behavior of the system, in anal- 
ogy to Ising spin systems, by computing the order param- 



FIG. 1. (Color online) Order parameter versus p for quenched 
(top) and annealed (bottom) random variables of Eq. 
with typical values of q indicated on the figure. The full 
(open) symbols are results of simulations with synchronous 
(asynchronous) update. The population size is N = 1000 and 
data are averaged over 100 realizations. 
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convictions evolve influenced only by the interaction with 
other randomly chosen agents. 

In Fig. [T]we exhibit results for the order parameter as a 
function of p for typical values of q, allowing to compare 
the cases with quenched (top) and annealed (bottom) 
disorder. One can see that the curves for synchronous 
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and asynchronous updates in the quenched case are al- 
most identical, indicating that the critical behavior is not 
modified by the update scheme when we consider frozen 
disorder. On the other hand, if we allow the disorder to 
fluctuate in time, the results for synchronous and asyn- 
chronous updates are different. 

We have verified numerically that, in all the analyzed 
cases of disorder and update scheme, the system under- 
goes a non-equilibrium phase transition for all values of 
q > 0. The transition separates an ordered phase, where 
one of the extreme opinions (+1 or —1) dominates, from a 
disordered one, where the three opinions coexist equally. 
A condition that was also obtained analytically for the 
synchronous annealed case (see the Appendix). 

In order to locate the critical points p c {q) numerically, 
we have performed simulations for different population 
sizes N . Thus, the transition points p c (q) are estimated, 
for each value of q, from the crossing of the Binder cumu- 
lant curves for the different sizes. In addition, a finite-size 
scaling analysis was performed, in order to obtain an esti- 
mate of the critical exponents (3, 7 and v. As an illustra- 
tion, we exhibit in Fig. [5]the behavior of the quantities of 
interest as well as the scaling plots for q — 0.5, quenched 
random couplings and synchronous updates. Our esti- 
mates for the critical exponents coincide with those for 
the original model (q = 1), i.e., we obtained /3 ~ 0.5, 
7 ~ 1.0 and \jv ~ 0.5. These exponents are robust: 
they are the same for all values of q, independently of 
the update scheme considered and of the kind of random 
variables (quenched or annealed). 

Taking into account the critical values p c (q) obtained 
from the simulations, we exhibit in Fig. [3] the phase dia- 
gram of the model in the plane p versus q. As discussed 
before, in the case of quenched variables the frontier is 




FIG. 2. (Color online) Order parameter, Binder cumulant 
and susceptibility for the PDF of Eq. Q with q — 0.5 and 
different population sizes N indicated on the figure (main 
plots). The corresponding scaling plots are shown in the re- 
spective insets. Data are for quenched random variables and 
synchronous update scheme. The best data collapse was ob- 
tained for p c ~ 0.167, ~ 0.5, 7 ~ 1.0 and 1/v ~ 0.5. 
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FIG. 3. (Color online) Phase diagram of the model defined 
by Eqs. (fl), (fSJl and ^ in the plane p versus q, separating 
the ordered (O) and disordered (D) phases. The symbols are 
the finite-size estimates of the critical points p c (q) obtained 
from the simulations, dashed lines are guides to the eye and 
the full line is the analytical result predicted by Eq. (|A2[) . 
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independent of the update scheme. On the other hand, 
for annealed variables the results are different. This is 
possibly due to the time fluctuation of the annealed vari- 
ables, which does not occur in the quenched case. The 
analytical prediction for the synchronous annealed case 
is presented in the Appendix. 

B. Model with distribution G2 

For the distribution Gi(Ci) of Eq. ((4]), a fraction 1 — q 
of the convictions d are —1 (instead of being null as 
in Sec. III. A). Now, the agents i present negative con- 
victions that contribute to a change in their opinions, 
together with the influence of a randomly chosen agent 

3- 

Differently from the case where the PDF of convictions 
is given by Eq. Q, analyzed in Sec. III. A, now we 




t 



FIG. 4. (Color online) Time evolution of the order parameter 
for N = 1000, p — 0.1 and typical values of q, labeling the 
curves from top to bottom. The curves are for: quenched 
asynchronous (top) and annealed synchronous (bottom). We 
can observe that a disordered state is reached when the value 
of q is decreased below a threshold. 



can observe that there is a threshold q c below which the 
system is always in a disordered state, for all values of p. 

The time evolution of the order parameter is illustrated 
in Fig. 2] for the quenched asynchronous and annealed 
synchronous cases. Similar evolution is observed for the 
other two combinations too. For sufficiently low q, none 
of the two extreme opinions dominates. Moreover, we 
verified that in such cases the fraction of each one of the 
three possible opinions is again 1/3 in average, indicating 
complete disorder. This result was also found theoreti- 
cally for the annealed synchronous case (see Appendix). 

In the cases where a transition occurs, a finite-size scal- 
ing analysis was perfomed as in Sec. III. A. The same 
mean-field exponents were obtained, independently of 
the update scheme considered and of the kind of ran- 
dom variables. The phase diagram for the different types 
of update and random variables is shown in Fig. [5] The 
analytical prediction for the synchronous annealed case, 
derived in the Appendix, is also included. In such case, 
Eq. (|A3|) predicts that for q < 3/4 no transition occurs, 
and the system is always in a disordered state. 
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FIG. 5. (Color online) Phase diagram of the model defined 
by Eqs. (TTJ, ((2| and ((4} in the plane p versus q, separating 
the ordered (O) and disordered (D) phases. The symbols are 
the finite-size estimates of the critical points p c (q) obtained 
from the simulations, the dashed lines are guides to the eye, 
and the full line is the analytical result given by Eq. (|A3|) . 
Notice that the transition is destroyed for values of q below a 
threshold, indicated schematically by arrows. 



IV. FINAL REMARKS 

In this work we have studied opinion dynamics through 
a model where agents interact by pairs in a fully- 
connected graph. Notice that it is related to the mean- 
field random-bond spin-1 Blume-Capel model [25[. In 
fact, the opinions have three different states (spin-1) and 
the agents interact through random couplings that can be 
positive/ferromagnetic (or negative/anti- ferromagnetic) 
with probability 1 — p (or p). Moreover, there is also a 



5 



self-interaction term, here interpreted as the conviction, 
while in the magnetic model it is associated to the crys- 
tal field. Then, models G\ and G2 are related to random 
diluted field and random antiferromagnetic impurities, 
respectively. The present kinetic model corresponds to 
zero temperature, but thermal fluctuations are substi- 
tuted by the randomicity of the couplings. We aimed to 
study the impact of this term on opinion formation and 
on critical behavior. Differently from other models of 
opinion formation, where the ordered state is marked by 
consensus, here the ordered state is characterized by the 
upraise of a dominating extreme opinion, which becomes 
consensual only in the limit in which interactions are all 
positive {p = 0). 

First, we have considered a PDF that aims to model 
populations where there is a fraction 1 — q of agents with- 
out self-convictions about their opinions, and thus they 
can be easily persuaded to change their opinions. Our 
results show that the critical fraction of negative inter- 
actions p c below which the population reaches partial 
agreement decreases smoothly for decreasing values of 
q, collapsing with p c = only at q = 0. This order- 
disorder transition is continuous and the critical expo- 
nents are universal and mean- field like, presenting the 
values (3 = 0.5, 7 = 1.0 and \ jv = 0.5 for all values of q. 

We have also considered a PDF for the convictions 
modeling societies where some agents have a tendency to 
change spontaneously their opinions. In this case, disor- 
dered states are favored, and the order-disorder boundary 
falls off rapidly to p c — for decreasing values of q. Thus, 
in opposition to the previous case, there are threshold 
values of q below which the system is always in the dis- 
ordered state. Despite this difference, the order-disorder 
transition is also continuous and the critical exponents 
are universal and mean-field like, as in the previous case. 
It is important to notice that the results depend quan- 
titatively (but no qualitatively) on the kind of update 
scheme used (synchronous or asynchronous) and on the 
nature (quenched or annealed) of the random variables 
considered, for the two studied PDFs. 



Appendix A 

Following the lines of Ref. [ll|, we computed critical 
values for the synchronous annealed case. We first ob- 
tained the matrix of transition probabilities whose ele- 
ments rriij furnish the probability that a state suffers 
the shift or change i — > j. Let us also define fx, fo and 
the stationary probabilities of each possible state. 

In the steady state, the fluxes into and out from a given 
state must balance. In particular, for the null state, one 
has mi.o + "J-1,0 = m o,i + "io,-i- Moreover, when the 
order parameter vanishes, it must be fx = f-\. In both 
cases considered below for the distribution of convictions, 
those two equalities imply fx = /_i = fo = 1/3 (disorder 
condition). This holds in particular at the critical point. 

Finally, let us define r(fc), with — 2 < k < 2, the prob- 



ability that the state shift per unit time is k, that is, 
r(k) = J2i m i,i+k- In the steady, the average shift must 
vanish, namely, 

2[r(2)-r(-2)]+r(l)-r(-l)=0. (Al) 



1. PDF Gi 

The elements of the transition matrix are 

mx,x = ff(l -p) + hfoq + fif-xP 

"M,o = ffqp + /i/o(l - q) + - p) 

mi,-! = A 2 (1 - q)P + - ff)(l -P) 

m ,x = io/i(l -p) + fof-xP 

m ,o = fo 
mo-x = /0/1P + /o/-i(l -p) 

m-1,1 = /i/-i(i - <?)(i - p) + f-xPi} - q) 
m-1,0 = fif-\q{i -p) + /o/-i(l - q) + f\qp 
"i-i-i = hf-ip + fof-iq + f-xO- -p)- 

The null average shift condition (|A1|) . together with 
the disorder condition, leads to 

•"° = W^y (A2) 

2. PDF G 2 

For this PDF, the transition matrix is 
"mx,x = /i 2 <?(l -p) + hfoq + fif-iqp 

mx,o = fx(qp + (1 - <7)(i - p)) + fxf-i{p + q - 2pq) 

mx,-x = A 2 (l - q)p + +/i/o(l -q) + - <?)(! -p) 

"io,i = /o/i(l -p) + fof-iP 

"10,0 = fo 

mo-x = fofxP + /o/-i(l -p) 

m_i,i = fxf-x(l - q)(l ~P) + /-i/o(l - q) + 7-iP(l - q) 
m-1,0 = /l/-l(p + q - 2pq) + f 2 _ x {{\ - p)(l -q) + qp) 

m-1,-1 = hf-ipq + faf-iq + /-i<?(i - p) • 

In this case Eq. (jAll) . toghether with the disorder 
condition, gives 

Pc = q - 3/4 . (A3) 

In contrast to the frontier defined by Eq. (|A2I) . which 
implies a critical value of p below which the system has a 
predominant opinion, Eq. (IA3|) implies that for q < 3/4 
the system can not order. 
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